Let’s create a bax plot showing the distribution of people’s order hour for each aisle
data("instacart")
box = instacart|>
group_by(aisle,order_hour_of_day)|>
mutate(n_order=n())|>
filter(n_order>1000)
plot_ly(box,x = ~aisle, y = ~order_hour_of_day,color=~aisle, type = "box",colors="viridis")